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ABSTRACT 


This thesis reports on three procedures and the 


associated numerical results for obtaining semiconductor 





majority carrier concentrations when subjected to a 


temperature sweep. The capability of predicting the 





exhaustion regime boundaries of a semiconductor is critical 








in understanding and exploiting the full potential of the 





modern integrated circuit. An efficient and reliable 


method is needed to accomplish this task. Silvaco 











International’s semiconductor simulation software was used 





to predict temperature dependent majority carrier 





concentration for a semiconductor cell. Comparisons with 








analytical and numerical MATLAB-based schemes were mad 
This was done for both Silicon and GaAs materials. 
Conditions of the simulations demonstrated effect known as 


Bandgap Narrowing. 
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EXCEUTIVE SUMMARY 


This thesis reports on three procedures and the asso- 


ciated numerical results for obtaining semiconductor 





majority carrier concentrations when subjected to a tempe- 
rature sweep. The capability of predicting the exhaustion 


regime boundaries of a semiconductor is critical in under- 





standing and exploiting the full potential of the modern 





integrated circuit. This includes, but not limited to, ICs 





designed for military hardware as well as specialized 





extrinsic semiconductors designed for space-based appli- 





cations. An efficient and reliable method is needed to 


accomplish the task of predicting exhaustion regime 





boundaries. Silvaco International’s semiconductor simula- 





tion software was used to predict temperature dependent 





majority carrier concentration for a semiconductor cell. 





Comparisons with analytical and numerical MATLAB-based 


schemes were made. Conditions of the simulations demon- 





strated th ffect known as Bandgap Narrowing. 








The goals of this project are twofold. The first goal 
is to design a program procedure to temperature sweep a 


uniform single-doped semiconductor with a two-dimensional 








Silvaco simulation software[1]. There is a significant 
challenge to overcome in meeting this goal. Specifically, 


this is due to a limitation in Silvaco software in that 








temperature sweeps cannot be performed within a single 


“run”. The second goal was to test the results obtained in 





Silvaco against one-dimensional analytical and numerical 
models implemented in MATLAB. 


Figure ES contains six plots of Jlogio(pPo) versus 


1000xX(1/T) for various conditions (with and without band- 


XV 


gap narrowing (BGN) and from various sources (numerical, 
analytic, and Silvaco) for the case of doped Silicon. The 
model for band-gap narrowing [6] used here is: 

ar’ 


E (T)=E,,-—— . 
1) =E,, (T+B) 





With and Without Bandgap Narrowing Plot 


19.0 














——1 Silvaco - w/o BGN 
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Figure ES.1 Composite Log Plot of Three Regimes of 
Silicon. 


The number codes identify the various curves as per 


description in the legend of the figure. 





In Table 1 


XV1 


recommended values for silicon including BGN parameters @ 





and B are provided. Note, as expected, curves from the 


numerical and analytic methods are indistinguishable, for 





example, MATLAB-based analytic curve #4(w/BGN) overlays 


MATLAB-based numerical curve #6(w/BGN). Furthermore, as 





xpected the impact of BGN is apparent in the curves at the 


higher temperatures. The difference between Silvaco pre- 











dictions and the 1-dimensional modeling efforts are slight 


but noted for future investigation. Similar analysis was 








completed for GaAs, with slight variations due to the 
different material structure and properties. The results 


for GaAs curves were Similar. 






































Parameter /| Value 
Ea, 0.045 eV 
Nz 10°’ #/cm 
Ego 1.11 eV (Silicon) 
cn 4 
BGN @ 0.00047 
BGN £B 636.0 
Table ES.1: Silicon Parameters used in Simulations 











TonyPlot was used to create a slideshow movie, to aid 





in visualization of the process. The structure files from 














the DeckBuild runs were saved and TonyPlot was used to 


display the hole/electron concentration of the lattice at 





temperature steps. This data was then exported into 





Microsoft Office PowerPoint. 
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i. INTRODUCTION 


A. CHAPTER INTRODUCTION 


This first chapter will provide an introduction to the 


work accomplished in this thesis, a literature review, and 





a structural overview of the thesis. 


B. PURPOSE 


The purpose of this research is to compare the 
performance of commercial semiconductor two-dimensional 


simulation tools with one-dimensional analytical and 








numerical models for freeze-out and exhaustion under a 


variety of nontrivial conditions. 


C. INTRODUCTION TO GENERAL TOPIC 








Integrated circuits (ICs) are specified to operate 


between designated temperature limits. The circuit designer 





selects the doping level or levels and typically assumes 





that the dopants are approximately 100% ionized, i.e., 





exhaustion and temperature are not too high. There can be a 
Significant impact on a plethora of device parameters such 


as depletions widths and/or FET threshold voltages if the 





assumption is violated. In the domain where the tempera 








ture is too low, the percentage ionization of dopant or 








dopants will be significantly less than 100%. The value 
for the majority carrier concentration is depressed 


significantly below the design value. On the other hand, if 





the temperature is too high, the thermal generation effect 
causes the majority carrier concentration to become 


excessively higher than the design value in what is called 
al 





the intrinsic temperature regime. Th xhaustion regime 


lies between this two extremes, intrinsic and freeze-out. 





It is well known that for a multiple impurity dopant 








process this "simple" thr regime description can be 





inadequate. What is important to the designer is the plot 





of the majority carrier versus the temperature or, what is 
more commonly done, a plot of majority carrier concen- 


tration versus the reciprocal of temperature. 


D. GOALS 


The goals of this thesis are twofold. The first goal 





was to design a program procedure to temperature sweep a 





uniform single-doped semiconductor using a two-dimensional 





Silvaco simulation software[1]. There is a significant 
challenge to overcome in meeting this goal. Specifically, 


there is a limitation in Silvaco software in that temp- 





erature sweeps cannot be performed within a single “run”. 








The second goal was to test the results obtained in Silvaco 








against one-dimensional analytical and numerical models 


implemented in MATLAB. 


E. THESIS OVERVIEW AND LITERATURE REVIEW 











This research begins in Chapter with a review of 








the analysis of the different portions of the project. The 




















analysis in Chapter serves as the background for the 




















MATLAB numerical and analytic methods. Chapter is an 
introduction and description for the Silvaco-based 


computational tools used to evaluate and compare with the 














MATLAB-based methods briefly covered in Chapter . The 


procedure for modeling of semiconductors in Silvaco and 





other technical aspects will be discussed. Chapter IV is 
2 


the presentation and discussion of the results obtained 





from simulations and testing. This data will be used to 








crosscheck the three approaches to obtaining temperature 








dependence of majority carrier concentration in singly 
doped semiconductors. Chapter V is a submittal of 
conclusions on the study completed and recommendations for 


further work. 





Several topics, which are discussed in the appendices, 
have been included for completeness. The corresponding 
information is presented this way, in order to avoid 
cluttering and overcrowding of the main text. Appendix A 


provides a brief introduction to general semiconductor 





principles. Appendix B provides an operational guide to 


reproducing the results in Silvaco, along with an 








abbreviated version of the Silvaco DeckBuild source code 


created for this work. Appendix C includes a pre-print of 














the paper presented at the 45° IEEE International Midwest 








Symposium on Circuits and Systems in Tulsa, OK, on Aug 4-7, 


2002. 


Currently, commercial software packages, such as 





available through Silvaco International [1], are well 


designed to solve the electron/hole transport problem. 





This type of calculation is usually required to predict the 








device I-V (current-voltage) characteristics. Surprisingly, 





using the same package to obtain a temperature dependent 








plot for majority carrier concentration for a uniform 





semiconductor requires a somewhat complicated procedure 


[2]. This thesis presents an efficient and novel way of 








obtaining this curve from the Silvaco International 








software and compares the results with a proposed one- 





dimensional single-eguation analytic model and a one- 





dimensional numerical model that predicts the temperature 


dependence for majority concentration in all regimes [3, 





4]. To-date, most "analytic" methods for determining 
dominant features in such plots make use of multiple ad-hoc 


arguments, which taken one at a time, applies in only two 











of the thr regimes mentioned above [5, 6, 7]. Numerical 


methods based on one-dimensional analysis provide 





significant flexibility in terms of making predictions when 





there are multiple dopants and when taking into account 








second order effects [4]. It can be argued that numerical 


methods are generally going to be applicable over a wider 





range of problems than what can be solved with exact 





analysis. Two and three-dimensional numerical finite 





element methods provided in CAD based commercial packages 
have been widely employed by the engineers in the 


semiconductor processing industry. The one-dimensional 





results, both numerical and analytical, are compared with 


predictions from such a package of tools from Silvaco 





International [1]. The numerical one-dimensional method is 





based on solving for zeros in the charge neutrality 


condition [5,6,7] iteratively for selected temperatures. 





The numerical algorithm involves a root finding scheme 








known as the Method of Interval Bisection [8]. The one- 





dimensional MATLAB-based numerical method can be applied on 





conditions not specifically studied in this thesis such as 





multiple dopants. Details on the one-dimensional numerical 





algorithm are not covered in this thesis and interested 





readers are encouraged to s References [3] and [4]. The 


analytic expression for the majority carrier temperature 


dependence is based on the well-known solutions to a cubic 


equation [9]. 


Band gap narrowing [6] is a well-known effect in which 


th ffective band gap of the semiconductor decreases with 





increasing temperature. Cases with and without the bandgap 





narrowing effect were simulated to exhibit the flexibility 





and consistency of the various approaches taken. Many of 








the salient features in this thesis have recently been 
summarized in a short paper [10]. Background information, 
which appears in Appendix A, was referenced from these 


sources: [11-15]. 


F. BENEFITS OF STUDY 


With the proliferation of semiconductor-based app- 


lications in modern society, it has become critical to 





understand the effects of temperature on semiconductors 








properties. Understanding the nature of temperature 
sensitivity before production and implementation of a 


device can reduce cost and risk factors for a final 








product. For example, the issue is particularly evident in 
the case of space-based applications, due to the harshness 


of the environment. 


G. CHAPTER SUMMARY 


This Chapter described the general area of 








investigation, the goals of the thesis, and provided a 





literature review. In the next chapter, the mathematical 





formulization, which serves as a basis for one-dimensional 


analytical and one-dimensional numerical scheme for solving 


for temperature-dependent majority carrier concentration, 


will be addressed. 


II. TECHNICAL ANALYSIS 


A. CHAPTER INTRODUCTION 


This chapter provides the background analysis for 
predictions of temperature dependent majority carrier 
concentration for both numerical and analytic MATLAB-based 


approaches (see Figure 2-1). 
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Figure 2-1: Layout of Project. 


B. BACKGROUND ANALYSIS 


The following assumptions were made: 1) The semi- 


conductor is spatially uniform, i.e., no spatial variation 





of the doping and energy bandgap. 2) The semiconductor is 





in equilibrium. 3) “Degeneracy condition” [6], which is a 
byproduct of excessive doping, has not occurred. These 
assumptions were used throughout the simulation for both 


MATLAB and Silvaco software. 





The Fermi-Dirac occupation probability for electrons 





and holes respectively are: 





and 
I - 1 
FD) = 1 — FOE) = ——aoare (2.2) 





where FE is the energy level defined on electron energy band 





diagram (see Figure 2-2). The Fermi energy level, Er, is 


, 





where electron and hole occupation probability are both 
50%, k is Boltzman’s constant, and T is temperature in 
Kelvin. Figure 2-2 is a representation of the energy 
bandgap diagram of a semiconductor. On the right hand side 


are the levels for the conduction and valance levels, E, and 


E,. The total bandgap size is represented as Ey, and AEg 





and AE, are used to represent the difference between th 


donor (Ng) and acceptor (Nz) levels and the band limits. 




















Figure 2-2. Diagram of Bandgap Layout. 
8 








The thr dimensional electron and hole density-of- 


states [6] are given by respectively: 


m 


ENS 72 
Ce) = 2( . (E — E.)''? H(E - E,) (253) 





hi? 
and 


— BE) (2.4) 


Vv 


V2 (m,° 
ra 


3/2 
(E, — Ey’? HE 





where H(x) is the Heavy-side function defined as: 


1 x 2 0 
H(x) = (2.5) 
0 otherwise 





and Ec, Ey are conduction and valance band edges defined on 








Figure 2-2, m, m, are th ffective masses for electrons 


e 


and holes at those band edges, and f is short-hand notation 





Planck’s constant divided by 27. 


The density of electron concentration in electron 


energy space is then obtained from combining Egs. (2.1) and 
(2.2): 
dn 
AHS SS ee A) (2.6) 
GiB) = = 


and, similarly, for holes using (2.2) and (2.3): 


ad = 
SPo = ¢ (E) * FD(E) Oe 


BF) = 
$,(E) i 


which is the density of hole-concentration in electron 


energy space. 


The equilibrium concentrations of electrons and holes 


are compactly expressed as: 


ny(Ep) = | 6,(E) dE (2.8) 


and 


P(Es) = [ 6. dE (2.9) 


—oo 


where implicit dependence on Fermi energy level is apparent 
through substitutions of Eqs. (2.1) through (2.7), into 
Eq.(2.8). The transformations, set as x=(E-E.)/kT for 
evaluation of no(Er) and set as x=(E,-F)/kT for evaluation 


of po(Er), leads to: 














* 3/2 o 1/2 
n(E,) = a8 ieee Va = f = ax (2.10) 
h 2? To ee 
and 
P(E.) = v2 (ime) va | 2 — (2.11) 
om Mu he 2 Vn 41 ier? 
where 
P= (Eo Bee (2.12) 
and 
y = (8, -— £,)/ kT. (2.13) 


The expressions defined by (2.10) and (2.11) can be 
simplified using a standard approximation [6] as shown in 


the next paragraph. 


The Fermi integral is given by [6]: 


Veg 
Fij,ly) = 7 (2A) 
+e 


a 


‘ 


For y < -3: 


Boy) ae [J fete -a| =e’, (2.15) 


This approximation applied in the context of evaluation of 





carrier concentrations, is known as the Boltzman’s Approx- 








imation. In physical terms, the applicability of the 





approximation depends on the Fermi level being located in 
the bandgap and not within 3kT of the conduction and 


valance band edges. 


After combining (2.10) thru (2.15) it follows that: 





(Gye Wee rey (25k6) 





and 
DP, (E,) = Nery. (2 pa) 
wher ffectiv lectron and hole density of states are 
respectively: 
2nm,kT) 
70 
N, = (=e) (2.18) 
and 
2am, kT) 
a 
N, = (==) (2.19) 


and Y, and WY, are degeneracy factors described in the next 
paragraph. 

The electron and hole concentration degeneracy factors 
[3] are: 


Let eS (25503 


for which y,=l for y<-3 and: 


| ee ee eo (2521) 
11 


fom white Yo ltor-y<=3. 


If the Boltzman’s Approximation can be applied, the 








usual case if the doping is not too heavy, then yY,~1 and/or 





X,~i. The standard Boltzman’s approximation form for the 
equilibrium carrier concentrations is recovered. Speci- 
fically 

nm 2 her re (222) 
and 

Py = Neer (2.23) 


Equations (2.22) and (2.23) can be taken as exact under 


intrinsic conditions, i.e., no impurity dopant employed. 


Two rules are needed to define the intrinsic 





condition. These rules predict the Fermi level under 








intrinsic conditions, Er;, and the extrinsic carrier 


concentration, n;. Under intrinsic conditions, electron and 





holes are created in pairs and therefore n,=n,=p,. The 


application using (2.22) and (2.23) leads to: 
n, = ./n,p, = .N.N,e 0/2" (2720) 


and after setting Er to Er; and solving, one determines: 


N Ee +k. 
E,, = kKT1n|—| + Be am (2329) 
N 2 








Cc 


Equations (2.24) and (2.25) can be expressed in term 


of the bandgap, F,, which is represented on Figure 2-2 as 





the energy difference between the conduction and valence 





band. If the valance band is used as a referenc nergy 








level, i.e., E,=0 as represented on Figure 2-2, E, is then: 
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Be Se ee (2.26) 


and it follows [6] that Eqs. (2.22) and (2.23) can be 


expressed in terms of n; and Er; as: 


pam (2527) 
and 

n, ea Se i ae (2.28) 
Equations. (2.27) and (2.28) are the equilibrium carrier 


concen-trations. 


Substitutional impurities can exist in one of two 








states. Specifically, ither the neutral state or an ion- 








ized state can exist. The basic rules that govern donor 





and acceptor ionizations [6] are: 











n° Ep ~ Eq 

a = ga } (2.29) 
Ny 

and 
nS 1 [== 
a it m), (2230) 
N, J. 


where subscript ‘d’ is for donor and ‘a’ is for acceptor, 
Na? and N° refer to neutral donor and acceptor impurity 

concentrations respectively, Nj and N, refer to ionized 
donor and acceptor concentrations, and gg and gz, are occu- 


pation indices. There ar lectron spin-based arguments 





[6,12] that indicate in most semiconductors g,=4 and g,=2. 


The fraction of ionization can be established via several 





observations. First, the sum of ionized parts, plus neu- 

















tral must be added to the total impurity used, specifically 


for donors: 
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N, = Ny + Ni (2231) 
and for acceptors: 


N, = NO +N. (32) 


a 


After dividing both sides of Eq. (2.31) by Ng it is 
found that: 





nt re n° 
[=| = aie < (23:33) 
N N 


where Ng and Nz are identified on Figure 2-2. 


Similarly, after dividing both sides of Eg. (2.32) by 
Nz, it is found that: 





~ 0 
[| Sg (2.34) 


Therefore, substitution of Eqs. (2.31 and 2.32) into Eqs. 
(2.33 and 2.34) leads to: 














N* 1 
dif —> (2503 54) 
N FE 'd 
i eee ( ED 
Fae 
and 
N_ 1 
FS = Ext) (2.36) 
i. 1+ ( ae 
Ge 


which provides the guidelines for Fermi energy level 





dependent degrees of ionization. 





Donor energy levels are typically cited using the 


conduction band as the reference. For similar reasons, the 





acceptor energy levels are typically cited using the 


valence band as a reference [3]: 





Ee B= AE, (2.37) 





and 


E, = AE, (2.38) 
which is consistent with Figure 2-2. 


Assuming we have Mg, mono-valent donors, and M:, mono- 
valent acceptors, the conditions of local charge neutrality 
exist. This is based on assumptions #1 and #2. For the 


condition: 


Ma a 


PU; B,) = py = ny (NG) S* (wz), =0 (2.39) 


1=1 j= 


where there is an implicit dependence on both Er and the 


temperature within each term of Eq. (2-22). The 





dependencies are obtained from Eqs. (2.35 to 2.38). For a 











specified temperature, Eq. (2.39) results in a non-linear 





Epi Ep) /kT 





expression in Z, where Z=e' At a specified Z, the 


Carrier concentrations can be predicted from: 











n, 
1 = (2.40) 
Z 
and 
oS ne (2.41) 
which is an equivalent revised version of Eqs. (2.27) and 
(2.28). 


Making use of the Fermi probability distributions to 








predict ionization levels of the donors and acceptors, 





specifically this requires the substitution of Eqs. (2.35, 








2.36, 2.40, and 2.41) into Eq. (2.39) and dividing by the 





intrinsic concentration which leads to a condition on 


P(Z,T): 
< Na ; N 
PZ)=0=Z+)>) 2 gt yy , (2.42) 
1+ K,2Z Pah 


The over-bar notation indicates division by n;. The tempera- 
ture dependent ionization factors: 


(Er: —Ea, ) / KT 


d = Gale : (2.43) 


x 
Ill 


and 


(Ea, -Ey;) / kT 


K, = 9,¢ (2.44) 


serves as a gauge for degrees of ionization of the associ- 





ated impurity. In summary, Eqs. (2.42 to 2.44) serve as a 





mathematical basis for both the analytic and numerical 





[3,4] MATLAB-based approaches. 


Cs ANALYTIC SOLUTION APPROACH 


This situation is described by the special case that 








M,=land M,=0. The summation subscript in Eq. (2.42) can 


be dropped to lighten the notation. Straightforward algebra 


leads to a cubic equation in @2Z: 


1 Nz a 
prie-(B|e-Leo, (2.45) 
a K 





The corresponding result for a single impurity N-type 
is exactly the same form, producing a cubic in Y=1/Z with 
revised coefficients obtained by letting the ‘a’ subscript 
(for acceptor) be replaced with the ‘d’ subscript (for 


donor). The solution can be defined in terms of 





coefficients for the reference cubic equation: 





x +ax +axt+a,=0. (2.46) 


To facilitate representation of a solution, the 








following intermediate parameterization of the problem is 


commonly taken [9]: 


Q = (3a, - af)/9 (2.47) 
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and 


R = (9aa, — 27a, — 2a) /54 (2.48) 

and 
D=Q +R (2.49) 
where, D, which is referred to as the “discriminant” for 





the cubic problem, will dictate the type of solutions 





possible. As per the fundamental theorem of algebra there 


will be three roots. For D negative, all roots are real 





while for D positive only one root is real while the other 








two are complex conjugates. If D is zero there will be 
repeated root. For the problem being evaluated it can be 
shown using symbolic mathematical methods that Q and D will 


always be negative. Furthermore, it can then be shown that 





of the three roots, only one is positive and, therefore, 





physically acceptable. That root will be predicted by the 
following formula [9]: 


Z = xX = 2j/-O cos (0 / 3) - 54 (23.0) 


where @ = cos” (r/J-0" ] once Z is determined, the majority 


carrier concentration, po, can be predicted from Eq. (2.41). 


D. CHAPTER SUMMARY 


This Chapter has established the mathematical basis 


for both MATLAB-based analytical and numerical predictor 








algorithms for temperature dependence of the majority 








carrier (see Figure 2-1). In this Chapter the main steps 





leading to a single equation analytic solution are 








described. In the next Chapter, the third approach, based 
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on the commercially available Silvaco International Soft- 


ware, is briefly explained. 
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III. SIMULATION 


A. CHAPTER INTRODUCTION 


For this project, a suite of commercially available 
simulation tools [1] intended for semiconductor testing was 


used. The primary objective of the thesis was to generate 








Silvaco predictions for the temperature dependent majority 
carrier. A comparison will be made with the MATLAB-based 
analytical and numerical methods discussed in the previous 


chapter (see Figure 2-1). The data from all three 





approaches were then compared graphically using Microsoft 





Excel. The following sections will discuss the processes 





developed in the simulation to generate the data for this 





comparison testing. 


B. SILVACO SIMULATION TOOLS 


The simulation in Silvaco was designed using the 





following modules from the suite: ATLAS, DeckBuild, 








TonyPlot and DevEdit. The following outline (Figure 3-1) 





illustrates the interconnecting relationships between the 





individual modules. 


19 













Structure 
File 





ATHENA 


Process 
Simulation. 


Command 
File 


Runtime 
Output 


Solution 
Files 











Figure 3-1. Silvaco 








Interaction Flowchart 


Ll 





(From 





ATLAS is a 2D/3D simulator for semiconductor devices. 


ATLAS can provide data and insights into the internal 


physical mechanisms of a device based on predicted 





electrical behavior. 











alone tool 


Fabrication environment. 


It can either be used as a stand- 





or aS a core unit for the Silvaco Virtual Wafer 


This module was used in 





connection with the other modules to deliver powerful and 


accurate data on the behavior of semiconductors 





[i]: 


DeckBuild is an interactive graphical user interface 


to provide a user-friendly runtime environment for 


integration of the different aspects of the Silvaco 


Software suite. 


creation and control. 


A control window is provided for file 


Many of the features are automated 





to allow for accurate simulation in a simple to use 








environment. 
can also be referenced in 


[elias 











generate a mesh for the structure designed in 


TonyPlot is the stand-alone program, 


which 





DeckBuild, to display the results 


DevEdit is a device editor which can be used to 


DeckBuild and 
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ATLAS. A limitation of device simulators prior to DevEdit 


was inadequate or poor structure meshes. DevEdit’s usage 





was integrated into DeckBuild to allow for a more complete 


and accurate solution [1]. 


Cc. DEVELOPMENT OF THE SILVACO SIMULATION PROCESS 





In order to develop the overall simulation, the first 





step was to understand the full scope of the problem and to 


determine the best features of Silvaco to find a solution. 





There were many automated features and defaults settings 


needed to be addressed to ensure that both sets of data 





were based on the same fundamentals. The DeckBuild coding 
in Silvaco followed the standard format of the program as 


illustrated in Figure 3-2: 





Group Statements 


1. Structure Specification ———— MESH 
REGION 
ELECTRODE 
DOPING 


—_—_————— MATERIAL 
MODELS 
CONTACTS 
INTERFACE 


Zr Mate ere Modells Nope crstlcd Eaton 


3. Numerical Models Specification —— METHOD 


4. Solution Specification LOG 
SOLVE 

LOAD 

SAVE 
De ReSsuLes  Anadlivicais ARNE 
TONYPLOT 














Figure 3-2. Input File Layout (After [1]). 
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In Figure 3-2, the left side, listed under Group, outlines 


how commands are clustered in Silvaco. The right side, 





listed under Statements, catalogs several of the commonly 


used commands and the order in which they are called. 











Several important values were set at the beginning of 














the DeckBuild code to ensure that default settings were not 





used. The specified parameters included: the dopant used, 





the e-alpha (Qa), e-beta (fB), and e-band gap (Ego) values, 


which were set before initializing the program. As pre- 





viously stated, the bandgap represents one of the most 








important parameters of a semiconductor and thus, it was 





set explicitly in DeckBuild coding. The equation for 
bandgap narrowing listed in the Silvaco ATLAS Manual [1] is 


as follows: 





(3.1) 


2, 2 
E(T) = E,(300) + a cas : | 


300+ B T+B 


where T is the temperature in degrees Kelvin and the @ and B 
parameters are semiconductor material dependent. £E,(300) is 
the energy gap at 300°K. An equivalent formulization used 
in the MATLAB-based program is in terms of E,(0) instead of 
Be(300)). [6 10) 





The structural dimensions of the device were then set 








and the contacts were placed with DevEdit. The device 


material was then set and the doping concentration was also 








set in DeckBuild. Choosing a model was the next challeng- 





ing aspect of research since certain DeckBuild coding 


options would decide the included physical effects. This 





was needed in order to obtain the most comprehensive 





solution to compare with the MATLAB-based approaches (see 


Figure 2-1). 
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For the model options in the Silvaco Simulations 





DeckBuild code, incomplete, conmob, and fermi were chosen, 


because they covered the widest comprehensive range of 








temperature-dependent effects. The incomplete portion was 





used to properly account for freeze-out when there were not 
sufficient amounts of thermal energy to assume 100% dopant 


Lonization. The conmob function was used to maintain 





concentration mobility as per the simple power rule for 


temperature dependence. The fermi command was used to 





account for cases of heavily doped materials. 


Along with the technical challenges of understanding 





the correct configuration for modeling a semiconductor, 














there were the difficulties in manipulating Silvaco due to 





inherent limitations of the program suite. A major 





obstacle related to the inability to program a temperature 





sweep in Silvaco within a single run. A method had to be 





developed to complete this process in an automated fashion. 








It was accomplished by replicating single-run code segments 








within DeckBuild. Each of the code segments was essen-— 





tially the same, except for the temperature, which was 





incremented consecutively (see Appendix B). The tempera- 





ture along with other pertinent data were recorded and 


logged after each step, which completed the ATLAS run. 








This process was repeated over a selected temperature 
range. Stored data could then be exported to a program, 


such as Microsoft Excel. This process allowed for 





significant flexibility in selecting the semiconductor 





parameters. It also allowed for an automation in the 





process which required very minimal user input. The 





overall process is illustrated in Figure 3-3. 
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Set Material Set Parameters Set Dimension 

































































(ie. Silicon, GaAs) ||(ie. E-alpha, Band gap) of device 
ATLAS RUN! 
Calculations based Probe of 
Set Temperature—» on inputs at [>| Temperature and 
specific temperature Carrier Conc. 
Calculation Saving 
of current |; of .log i. é rene 
and voltages and .str files 























| 


NEXT ATLAS RUN! 























Figure 3-3. Flowchart of Silvaco Process. 


D. DEVELOPMENT OF THE MOVIE 


TonyPlot was used to create a slideshow, to aid in 
visualization of the process. The structure files from the 


DeckBuild runs were saved and TonyPlot was used to display 





the hole/electron concentrations of the lattice over the 
temperature steps. These data outputs were then exported 
into other commercial software packages that allowed for 


graphic file animation, such as Microsoft Office Power- 





Point. Two samples of the movie are provided in Figure 3-4 


and 3-5. The arrow over the curves on the right side iden- 





tifies the temperature. The upper half shows a selected 
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movie slide when the temperature was approximately 10 K 


while the lower half was taken at a much higher temperature 





of 900 K. Note that colors for the left side are legend 





encoded for the concentrations of the majority carrier at 





the corresponding temperature. More details on the 








procedure for creating a movie are provided in Appendix B. 


E. CHAPTER SUMMARY 


This Chapter outlined and explained the technical 


aspects of utilizing commercially available Silvaco 








International software for the objectives of this thesis. 








A more detailed description of this procedure is provided 





in Appendix B. The next chapter will provide a review and 





comparison of the Silvaco based results with predictions 





generated from MATLAB-based analytic and numerical 





algorithm. 
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IV. RESULTS ANALYSIS AND COMPARISON 


A. CHAPTER INTRODUCTION 


This chapter provides a comparison of the results for 
various approaches obtained from examples taken of Silicon 
and GaAs. The different conditions and cases documented 
will be discussed. For this project, Silvaco Suite was 


used to develop a procedure to predict temperature 





dependence of majority carrier concentration. This was 





discussed in some detail in the previous chapter. The 


results from Silvaco are compared with MATLAB-based 





analytic and numerical predictions. The data from both 





simulations are then compared graphically in Microsoft 





Excel. The following sections will discuss the processes 








developed in the simulation to generate the data for this 





comparison testing. 


B. RESULTS FOR SILICON 


Figure 4-1 contains six plots of l10o0gi9 (po) versus 
1000*(1/T) for various conditions from the various sources 
(numerical, analytic, and Silvaco). The data was also 
grouped into the case with and without bandgap narrowing 
(BGN). The model for bandgap narrowing [6] used here is: 

ar* 


E(T) = E,. - —————_. 4.1 
Ae) 2 e+ B) ( ) 


Numbers codes for the various curves on Figure 4-1 appear 


with the description in the legend. Recommended values 


[2,15] for silicon, including BGN parameters @ and Bf, appear 





in Table 4-1. Note, as expected, curves from the MATLAB- 
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based numerical and analytic methods are indistinguishable. 





For example, MATLAB-based analytic curve #4(w/BGN) overlays 


MATLAB-based numerical curve #6(w/BGN). Furthermore, the 





impact of BGN is apparent in the curves at the higher 


temperatures. 





With and Without Bandgap Narrowing Plot 





—1 Silvaco - w/o BGN 


18.5 Silvaco - w/ BGN 





Analytic —- w/o BGN 


18.0 Analytic - w/ BGN 





Numerical - w/o BGN 












































17.5 Numerical - w/ BGN 
fe) 
Q, 3,4,5&6 
“ 17.0 / 
ad 
fo) 
fe] 
4 
16.5 
16.0 
15.5 
15.0 
0.50 2.50 4.50 6.50 8.50 10.50 
1000 X 1/T (K) 
Figure 4-1. Composite Log Plot of three Regimes of 


Silicon. 
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The differenc between Silvaco predictions and the 
one-dimensional modeling efforts are slight but noted for 


future investigation. 






































Parameter Value 

Ez 0.045 ev 

N, 10°’ #/cm? 

Ego 1.11 eV (silicon) 

Ja 4 

BGN @ 0.00047 

BGN B 636.0 

Table 4-1. Silicon Parameters. 
Cc. RESULTS FOR GALLIUM ARSENIDE 








The results for GaAs seen in Figure 4-2 were similar 


to results appearing in Figure 4-1. The parameters used 





for this test are provided in Table 4-2 [2, 15]: 









































Parameter Value 

Ea 0.15 eV 

N, 10°* #/cm* 

Ego 1.42 eV (GaAs) 
Ja 4 

BGN @ 0.00054 

BGN £B 204.0 

Table 4-2. GaAs Parameters. 
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With and Without Bandgap Narrowing Plot for GaAs 





— 1 Silvaco w/o BGN 





— 2 Silvaco w/ BGN 





—— 3 Analytical w/o BGN 





— 4 Analytical w/ BGN 


— 5 Numerical w/o BGN 





Logi (PO) 


— 6 Numerical w/ BGN 




















0.00 2.00 4.00 6.00 8.00 10.00 12.00 
1000 X 1/T (K) 











Figure 4-2. Composite Log Plot of three Regimes of 
GaAs. 


Again, all three regimes, i.e., freeze-out, exhaustion, and 
intrinsic, were evident on Figure 4-2. As before, in the 
exhaustion regime, the majority carrier was approximately 


equal to the doping level, i.e., 10% #/cm*. The curves 





indicated that the MATLAB-based approaches produced pre- 








dictions that were indistinguishable. It should also be 





noted that the level of agreement between the Silvaco-based 
prediction and the MATLAB-based approaches was better in 


the GaAs case (see Figure 4-2) versus the silicon case (see 
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Figure 4-1). This is especially true in comparing the 


level of agreement in the freeze-out regimes of both 





Figures 4-1 and 4-2. As it turns out, the testing pro- 


cedure for the two materials, Si and GaAs, was essentially 








the same with one important exception. For the Silvaco- 
based simulation of Silicon, the input specification for 


the acceptor energy level (Fz) was not done directly. 











Instead, the chemical impurity dopant Boron was specified 





in the DeckBuild code. In the case of the GaAs test, this 











feature was overridden by a direct specification of the 








acceptor energy level (Fz) in the DeckBuild code. The 











latter approach is clearly aligned with the MATLAB-based 

















methods. It is a reasonable conjecture that this accounts 
for the noticeable improvement in the level of agreement 


between the Silvaco and MATLAB-based methods. 


D. DISCUSSION OF RESULTS 





The level of agreement between results generated from 


the MATLAB approaches and from the Silvaco source code 





tended to confirm the accuracy of the models. The Silvaco 


and MATLAB-based approaches were in close agreement in 





predictions for the majority the carrier concentration over 


all three single impurity regimes of freeze-out, exhaus-— 





tion, and intrinsic. 


An important effect included in the testing was 
Bandgap Narrowing (BGN) [6]. The results were categorized 
for each material tested with BGN and without BGN. The 


differences observed for BGN proved to be fairly 





Significant at higher temperatures. In the exhaustion and 


freeze-out regimes, the impact of BGN was not noticeable. 
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E. CHAPTER SUMMARY 


Both the one-dimensional analytic model and numerical 





algorithm were essentially indistinguishable in terms of 








predictions. As expected, when BGN effects were included 
in the models, the impact was observed at the higher 


temperatures. The difference between predictions of the 





Silvaco-based two-dimensional modeling and the one- 








dimensional solutions were observed to be slight but worthy 
of future investigation. The next Chapter will summarize 
the main conclusions and will provide suggestions for 


future work. 
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Vv. CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY OF WORK 


This Chapter provides a summary of the thesis and 
provides some final thoughts on future areas of study and 


improvement. The principal completed goal of this research 





has been to design a convenient automated procedure that 





will predict, via Silvaco International semiconductor 





simulation software [1], temperature-dependent exhaustion 





and freeze-out effects in single-dopant, homojunction 





semiconductors. To gauge the level of success of this 
effort the results generated from the Silvaco simulation 


software have been compared with similar predictions from 





MATLAB-based analytical and numerical algorithms. This was 





done for both silicon and gallium arsenide semiconductors. 





Chapter I provided an introduction to the objectives 





and benefits of this research, along with a preview of what 














topics would be addressed. Chapter presented an in- 





depth mathematical analysis, which serves as a basis for 





MATLAB-based analytic and numerical approach. This chapter 


also detailed the process used to determine the analytical 























solution. Chapter was focused on imparting an overview 





of Silvaco and how it was implemented to provide the neces 





sary results. Chapter IV was used to present the data 
acquired from the different tests and to make a comparative 


study based on the data. A sample selection of slides from 





the Silvaco movie prepared during the thesis is also 





presented in Chapter IV. 
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B. CONCLUSION 


This thesis was designed to be part of a larger study 


in creating very robust MATLAB-based analytic and numerical 





simulations to provide an accurate portrait of temperature 
effects occurring in a semiconductor. This would allow for 
a versatile and reliable method in predicting the regimes 


of a semiconductor. 


After reviewing the data, it can be concluded that the 
MATLAB-based analytical and numerical models compared very 


favorable to the results obtained in Silvaco. All three 








methods accurately predict the semiconductor behavior in 


the three regions from freeze-out to the intrinsic. The 





MATLAB-based analytic and numerical coding allowed for a 





quick and robust method to generate the majority carrier 
concentration data of a semiconductor under a varying 


degree of changes in material and doping level. Currently, 





the computation cost of results obtained via MATLAB 











required much less time than the Silvaco method by roughly 


a factor of 50. 


Cc. RECOMMENDATIONS 


There has been significant progress made in this 





research, but there still exist much additional work to be 
accomplished in this area. The future of analytical 
modeling in MATLAB needs to be refined and expanded to 
allow for a greater range of effects and conditions. For 


example, the range of cases that can be solved analytically 





could be widened to include multiple impurities. This 





capability already exists in the MATLAB-based numerical 





algorithm. At this present time, the numerical algorithm 
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has not satisfactorily been tested for high levels of 





doping when the Boltzman’s approximation fails. 





There still exist many conditions and parameters in 





Silvaco which have not been accounted for which could be 
implemented in a faster MATLAB model. A more efficient 
method of handling temperature ranges in Silvaco is also 


something that needs to be improved upon. 


Finally, utilizing Silvaco to create and test a three- 











dimensional model might provide more conclusive data to 
continue to validate the results obtained for the 


analytical model. This would provide a significant insight 








into the limits associated with the one-dimensional models 


being tested in MATLAB. 
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APPENDIX A. INTRODUCTION TO GENERAL SEMICONDUCTOR 
PRINCIPLES (AFTER [11, 12]) 


A. QUANTUM MECHANICS 


To fully understand the basis of this thesis, a basic 


knowledge of quantum effects is required to understand why 





particles behave as they do at the sub-atomic and atomic 








dimensional scale. Before the twentieth century, scien- 











tists had a firm belief that the physical world was divided 


into two worlds, the world of particles and the world of 





waves. This belief was based on the classical system 





theorized by Sir Isaac Newton. 





In the 20° Century, this particle/wave paradigm was 
challenged and proven to be incomplete, based on the early 


work of physicist such as Neil Bohr. Werner Heisenberg, 





and an Austrian, Erwin Schrédinger postulated a new Quantum 
Theory based on the aforementioned work, and thus a valid 
and proven method was constructed to provide a means of 


describing sub-atomic particles. 


Quantum mechanics theorizes that the electronic struc- 


ture of an atom is based on quantum conditions. These 








conditions are postulated in the theory that electrons in 





an atom are restricted to certain energy level. This 





“quantization” defines certain allowable transitions invol- 
ving absorption and emission of energy by the electrons 


ise a 


Neils Bohr determined that electrons in orbits moved 





with well defined momentum. Based on this work and through 





absorption and emission, an electron may transition to the 
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next smaller or larger energy level based on the following 














rules: 

1) Electrons exist in certain stable, circular 
orbits about the nucleus. 

2) The electron may shift to an orbit of higher 
or lower energy, thereby losing or gaining 
energy equal to the difference in the energy 
levels (by absorption or emission of a photon 
of energy hv, see Figure A-1). 

3) The angular momentum Pg of the electron in an 
orbit is always an integral multiple of 
Planck’s constant divided by 27m [11]. 

E2 
aS. 
WEL 
Figure A-l. Emission of Photon. 
B. SEMICONDUCTOR MATERIAL PROPERTIES 


A semiconductor is able to conduct electricity at room 





temperature more readily than an insulator, but less easily 











than a metal. Electrical conductivity, which is the 
ability to conduct electrical current under the application 
of a voltage, has one of the widest ranges of values of any 
physical property of matter. Such metals as copper, 


Silver, 





and aluminum are excellent conductors, but such 
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insulators as diamond and glass are very poor conductors. 
At low temperatures, pure semiconductors behave like insu- 
lators. Under higher temperatures or light or with the 
addition of impurities, however, the conductivity of 


semiconductors can be increased dramatically, reaching 








levels that may approach those of metals [12]. 


The common semiconductors include chemical elements 
and compounds such as silicon, germanium, gallium arsenide, 
and lead telluride. They increase in conductivity with 
temperature, light, or impurities arising from an increase 


in the number of conduction electrons, which are the 





carriers of the electrical current. In a pure, or intrin- 





sic, semiconductor such as silicon, the valenc lectrons, 





or outer electrons, of an atom are paired and shared bet-— 


ween atoms to make a covalent bond that holds the crystal 





together. These valenc lectrons are not free to carry 





electrical current. To produce conduction electrons, 





temperature or light is used to excite the valenc lec 





trons out of their bonds, leaving them free to conduct 











current. Deficiencies, or “holes," are left behind that 
contribute to the flow of electricity. (These holes ar 
said to be carriers of positive electricity.) This is the 





physical origin of the increase in the electrical conduc-— 








tivity of semiconductors with temperature. The energy 
required to excite th lectron and hole is called the 
energy gap. 


For a semiconductor, this energy gap is termed the 





bandgap. It is a forbidden region in which an electron 











cannot exist, as it is the “space” between the quantized 


states. The upper region is called the conduction band and 
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the lower is labeled the valence band. This separation 


between th nergy of the two bands called the bandgap is 








labeled as E,. It is the most important parameter when 


dealing with semiconductors since it is the basis for any 





utilization of a given material. These bandgaps define the 
relative ease to “free” electron through excitation from 
thermal, optical, or electrical energy. Semiconductors 
have lower bandgaps (typically 0.5-4 eV) while insulators 
have higher bandgaps (5-10 eV) of transition. 


Semiconductor devices have many varied applications in 





electrical engineering. Recent engineering developments 











have yielded small semiconductor chips containing millions 


and in the near future, hundreds of millions, of transis-— 





tors. These chips have made possible great miniaturization 





of electronic devices. More efficient use of such chips 





has been developed through what is called complementary 





metal-oxide semiconductor circuitry, or CMOS, which con- 
sists of pairs of p- and n-channel transistors controlled 


by a single circuit. 





Cc. CONDUCTION 


Certain critical terms and concepts must be introduced 








to better understand the fundamental of semiconductors. 
Two parameters in the handling and usage of semiconductors 


are temperature and doping. These two factors are extreme- 





ly important in determining how effect a semiconductor is 





and how it can most effectively be used in any system 





design. 


Doping is a method to produce free carriers of elec-— 





tricity in a semiconductor by adding impurities to, or to 
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"dope," the semiconductor. The difference in the number of 
valence electrons between the doping material, or dopant 
(either donors or acceptors of electrons), and host gives 


rise to negative (n-type) or positive (p-type) carriers of 





electricity. This concept is illustrated in the accom-— 
panying diagram of a doped silicon (Si) crystal. Each 

Silicon atom has four valence electrons (represented by 
dots in Figure A-2); two are required to form a covalent 


bond. In n-type silicon, atoms such as phosphorus (P) with 








five valence electrons replace some silicon and provide 








extra negative electrons. In p-type silicon, atoms with 





three valence electrons such as aluminum (Al) lead to a 








deficiency of electrons, or to holes, which act as positive 


electrons. The extra electrons or holes can conduct 





electricity. 





7 2 be 7 be 7 i 7 be 


a 
“Sr Sr ° Po Sb Sr electron 


“$0 $0 “Sr Al Sr 


7 7 be 7 7 be 7 hole 














Figure A-2. Si Material Doped with P (From [13]). 


When p-type and n-type semiconductor regions are 
adjacent to each other, they form a semiconductor diode, 
and the region of contact is called a p-n junction. (A 
diode is a two-terminal device that has a high resistance 
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to electric current in one direction but a low resistance 








in the other direction) The conductance properties of the 





p-n junction depend on the direction of the voltage, which 








can, in turn, be used to control the electrical nature of 








the device. Series of such junctions are used to make 








transistors and other semiconductor devices such as solar 


cells, p-n junction lasers, rectifiers, and many others. 


Temperature is also very important and as it is raised 
from O°K, electrons are excited from the thermal energy and 


they move across the bandgap from the valence band to the 





conduction band. The resulting material is one with some 
unoccupied states in the valence band and this vacancy 


creates a “hole”, or an empty state. This creates an 








electron-hole pair (EHP) and as these “freed” electrons 
move through the material, this motion induces an electric 


field and current flows. 


There are two major classes of semiconductors, intrin- 
sic and extrinsic. An intrinsic semiconductor is one in 


which the electrical properties of device is determined by 








the host material, such as with bulk silicon. In an 


extrinsic semiconductor, the chemical impurities used to 








dope the host material are the determining factor for the 


electrical properties of the material. These dopants can 








drastically modify the electrical properties of a device as 











with gallium arsenide. 
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Figure A-3. Concentration vs. Temperature. 





As illustrated in Figure A-3, when the temperature 
increases, the intrinsic properties of the material over- 


ride any level of doping applied to the material. This 





region varies based on material and doping level of im- 





purities. At low temperatures, the energy from thermal 
effects is not high enough to fully activate all the donor 
and acceptor impurities. This region is known as the 
freeze-out region and this topic will be addressed later in 


this report. 





If temperature increases above 0 K, excitation of 


electrons is induced and EHPs are formed as thermally 





excited electrons cross the bandgap to the conduction band. 


The positions of the free electron and hole in the crystal 
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are spread out over several lattice spacing and can be 





considered quantum mechanically by probability distri- 





butions. At a given temperature, there is a certain amount 





of EHPs. In a steady state condition the generation rate 
(G;) of EHPsS in an intrinsic material must be equal to the 


rate of recombination of electrons and holes [11]. 


As temperature increases, the probability distribution 





of an electron in a range of allowed energy at equilibrium 





follows Fermi-Dirac statistics, because they obey the Pauli 











exclusion principal. The Fermi-Dirac distribution function 
is equal to: 


1 
JE) = er (he 1) 


where E£, is called the Fermi Level. For intrinsic ma- 
terials, the probability an electron is at the Fermi level 
is very near %. When T=0K, this function is a unit step 
and, as T increases, the function takes on a more curved 


appearance (see Figure A-4). 
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Figure A-4. Fermi-Dirac Distribution. 
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APPENDIX B. OPERATIONAL GUIDE TO REPRODUCING THE 
RESULTS IN SILVACO 


A. OPERATIONAL GUIDE TO REPRODUCING SILVACO DATA 
ACQUISITION CAPABILITIES 
Silvaco’s semiconductor software is a powerful and 

very capable tool in deciphering the electrical and 

material properties of a semiconductor. The vast array of 


built-in functions and add-on modules create a powerful 





tool, albeit a tool with a steep learning curve. This 














concise section will attempt to help create a standardize 








method to collect, organize, and extract the data from 


Silvaco for further analysis. 





The process begins with the DeckBuild graphical user 





interface and the creation of an input file. ATLAS is the 





module used in this research and the module is adequate and 





effective for all of the work that was accomplished. The 
physical parameters and structural dimensions were set for 


the input file; to be used globally throughout the model, 





in conjunction with a doping level and MODEL command 





setting. In the MODEL command, the lattice temperature is 








set for the ATLAS run. At this point, the program is ready 





to enter the recording stage of the process. 








The PROBE function was used extensively to gather the 


data required for proper analysis. This function was used 








to record the temperature, the hole concentration, and the 





electron concentration. An iterative stepping of minimal 











bias is applied to allow for the software to calculate the 





properties of the material. The STRUCTURE (.STR) file is 





saved and the probed values are logged into a LOG (.1log) 
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file. At this point, the program is ready to quit the 





ATLAS run at the given temperature and begin stepping on to 





the next temperature to begin the process anew. The new 











temperature is set in the MODEL command. The LOG file is 


then appended to record the new data from the updated 





temperature. The process repeats itself until the QUIT 
command is called to end the process. This raw data can 
then be exported into Microsoft Excel or MATLAB to be 


further categorized and analyzed. 


B. DECKBUILD ABBREVIATED SOURCE CODE 


go atlas 

set dopant = boron 

# Filneame 

set datadump = ThesisSl 
set e-alpha = 0.00047 


set e-beta = 636 


set e-band = 1.42 


0.045 


set Kae 





set Pae = 0.02 
mesh space.mult=1 


x.mesh loc = 0.00 spac = 0.10 





x.mesh loc = 1 spac = 0.10 
y.-mesh loc = 0.00 spac = 0.02 
y.-mesh loc = 1 spac = 0.02 


# placing anodes and cathodes 
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region number = 1 material = Si 


elec 
elec 








material reg 
EG300=S"e-ba 


trode name = anode top 
trode name = cathode bottom 


ion=1 egalpha=$"e-alpha" egbeta=$"e-beta" 








doping reg = 


model temp = 








nda" FAB =S Ww Fae" EDB =S "Dae" 





uniform $"dopant™" conc=lel7 





100 incomplete conmob fermi 





PROBE NAME = 
PROBE NAME = 


y = 0.2 





Temperature_(k) LAT.TEMP x = 0.05 y = 0.2 
Hole Concentration_(1/cm*3) P.CONC x 








PROBE NAME = 
0.05 y = 0.2 


sol 
sol 
sol 


lve init 
lve vanode 
lve vanode 





save outf = 
save outf = 
log outfile = S$"datadump".log 


solve 


go 





atlas 


Electron _Concentration_(1/cm*3) N.CONC 


= 0.1 

= 0.2 
S"datadump".str 
S"datadump"60.str 





vanode = 1 


mesh infile=S"datadump".str 





material reg 
EG300=$"e-band" FAB =S"Eae" EDB =S"Pae" 


models temp 


ion=1 egalpha=$"e-alpha" egbeta=$"e-beta" 











= 110 incomplete conmob fermi 





PROBE NAMF= 
PROBE NAMF= 
y=0.2 


Temperature_(k) LAT.TEMP x=0.05 y=0.2 


Hole Concentration_(1/cm*3) P.CONC x=0.05 








PROBE NAME = 
0.05 y = 0.2 


sol 





Electron _Concentration_(1/cm*3) N.CONC 





lve init 


save outf=S"datadump".str 


sol 
sol 


lve vanode 
lve vanode 





save outf = 


= 0.1 
= 0.2 
S"datadump"110.str 
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x 





log outfile=S"datadump".log append 
solve vanode=1 

go atlas 

mesh infile=S"datadump".str 


material region=1 egalpha=S"e-alpha" egbeta=S$"e-beta" 
EG300=S"e-band" EAB =S"Rae" EDB =S"Pae" 














models temp = 120 incomplete conmob fermi 





PROBE NAME= Temperature_(k) LAT.TEMP x=0.05 y=0.2 








PROBE NAME= Hole Concentration_(1/cm*3) P.CONC x=0.05 
y=0.2 
PROBE NAME Electron_Concentration_(1/cm*3) N.CONC x 








0.05 y = 0.2 


solve init 








save outf=S"datadump".str 
solve vanode = 0.1 

solve vanode = 0.2 

save outf = S"datadump"120.str 








log outfile=S"datadump".log append 
solve vanode=1 





go atlas 
mesh infile=S"datadump".str 


material region=1 egalpha=S"e-alpha" egbeta=$"e-beta" 
EG300=S"e-band" EAB =S"KRae" EDB =S$"Pae" 














models temp = 1010 incomplete conmob fermi 





PROBE NAME= Temperature_(k) LAT.TEMP x=0.05 y=0.2 





PROBE NAME= Hole Concentration_(1/cm*3) P.CONC x=0.05 
y=0.2 
PROBE NAME Electron _Concentration_(1/cm*3) N.CONC x 











0.05 y = 0.2 


solve init 
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save outf=S"datadump".str 

solve vanode = 0.1 

solve vanode = 0.2 

save outf = S$"datadump"1010.str 
log outfile=S"datadump".log append 
solve vanode=1 











go atlas 
mesh infile=S"datadump".str 


material region=1 egalpha=S"e-alpha" egbeta=S$"e-beta" 
EG300=S"e-band" EAB =S"Kae" EDB =S"Pae" 














models temp = 1020 incomplete conmob fermi 





PROBE NAME= Temperature_(k) LAT.TEMP x=0.05 y=0.2 








PROBE NAME= Hole Concentration_(1/cm*3) P.CONC x=0.05 
y=0.2 
PROBE NAME Electron_Concentration_(1/cm*3) N.CONC x 








0.05 y = 0.2 


solve init 

save outf=S"datadump".str 

solve vanode = 0.1 

solve vanode = 0.2 

save outf = S"datadump"1020.str 
log outfile=S"datadump".log append 
solve vanode=1 











go atlas 
mesh infile=S"datadump".str 


material region=1 egalpha=S"e-alpha" egbeta=$"e-beta" 
EG300=S"e-band" EAB =S"Rae" EDB =S$"Pae" 














models temp = 1030 incomplete conmob fermi 





PROBE NAME= Temperature_(k) LAT.TEMP x=0.05 y=0.2 





PROBE NAME= Hole Concentration_(1/cm*3) P.CONC x=0.05 
y=0.2 
PROBE NAME Electron _Concentration_(1/cm*3) N.CONC x 











0.05 y = 0.2 


solve init 
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save outf=S"datadump".str 

solve vanode = 0.1 

solve vanode = 0.2 

save outf = $"datadump"1030.str 
log outfile=S"datadump".log append 
solve vanode=1 











go atlas 
mesh infile=S"datadump".str 


material region=1 egalpha=S"e-alpha" egbeta=S$"e-beta" 
EG300=S"e-band" EAB =S"Rae" EDB =S"Pae" 














models temp = 1040 incomplete conmob fermi 





PROBE NAME= Temperature_(k) LAT.TEMP x=0.05 y=0.2 








PROBE NAME= Hole Concentration_(1/cm*3) P.CONC x=0.05 
y=0.2 
PROBE NAME Electron_Concentration_(1/cm*3) N.CONC x 








0.05 y = 0.2 


solve init 

save outf=S"datadump".str 

solve vanode = 0.1 

solve vanode = 0.2 

save outf = S"datadump"1040.str 
log outfile=S"datadump".log append 
solve vanode=1 











go atlas 
mesh infile=S"datadump".str 


material region=1 egalpha=S"e-alpha" egbeta=$"e-beta" 
EG300=S"e-band" EAB =S"KRae" EDB =S$"Pae" 














models temp = 1050 incomplete conmob fermi 





PROBE NAME= Temperature_(k) LAT.TEMP x=0.05 y=0.2 





PROBE NAME= Hole Concentration_(1/cm*3) P.CONC x=0.05 
y=0.2 
PROBE NAME Electron _Concentration_(1/cm*3) N.CONC x 











0.05 y = 0.2 


solve init 
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save outf=S"datadump".str 

solve vanode = 0.1 

solve vanode = 0.2 

save outf = $"datadump"1050.str 
log outfile=S"datadump".log append 
solve vanode=1 











Cc. PROCEDURAL GUIDE TO CREATING A MOVIE 


The ability to visualize the effects of temperature on 
a semiconductor cell is a very impressive method for illus- 
trating the concepts involved. For this purpose, Silvaco 
has the capability to save STRUCTURE files from the 


DeckBuild computations. These files are repositories of 





information for the semiconductor cell at given conditions. 





These files can be used to display the carrier concentra- 





tion, along with a host of other material and electrical 


properties. To access these files, Silvaco provides a 





stand-alone viewer called TONYPLOT. 





Multiple STRUCTURE files can be loaded into TONYPLOT, 


and the display drop down menu item can be used to display 





the carrier concentration of the files. All of these files 





must be selected because the programs sets the displayed 

















picture relative to the structure windows selected. While 
there is a effective movie animator built into Silvaco, the 
created movie is not exportable into other media format, 


thus each individual file must be screen captured into a 





standard graphical file type such as (.jpg) or (.bmp), for 
future manipulation. The files can then be organized in 
commercial moviemaker software and formatted for further 
animation. (Microsoft Office PowerPoint was used to animate 


the slides into a movie.) The process is labor intensive, 
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but very effective in creating an illustrative tutorial of 


the desired simulation effect. 
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APPENDIX C. PAPER PRESENTED AND ACCEPTED AT THE 
45°" IEEE INTERNATIONAL MIDWEST SYMPOSIUM ON 
CIRCUITS AND SYSTEMS IN TULSA, OK ON AUG 4-7, 2002 


This appendix contains the “Comparison of Analytic and 


Numerical Models with Commercially Available Simulation 








Tools for the Prediction of Semiconductor Freeze-out 


Exhaustion” paper from the 4-7 August 2002 45*® IEEE Midwest 





Symposium on Circuits and System at Oklahoma State 


University in Tulsa, Oklahoma. The paper was a condensed 





report on the work researched in this thesis. 
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Abstract - Currently, commercial 
software packages, such as 
available through Silvaco 
International[1], are well 
designed to solve the 


electron/hole transport problem. 


This type of calculation is 
usually required to predict the 
device IV characteristic. 
Surprisingly, using the same 
package, to obtain a temperature 
dependent plot for majority 
carrier concentration for a 


uniform semiconductor requires a 
somewhat complicated procedure 
[2]. Our paper will present an 
efficient novel way of obtaining 


this curve from the Silvaco 
International software and 
compare the results with a 


proposed one dimensional single- 
equation analytic model and a 
numerical model that predict the 
temperature dependence for 
Majority concentration in all 
regimes [3] [4]. 

1. Introduction 
Integrated Circuits 
specified to operat between 
designated temperatur limits. 


The circuit designer selects the 
doping level or levels and 


are 











56 





typically assumes that the 
dopants are approximately 100% 
ionized, i.e. exhaustion and the 
temperatur ar not too high. 
There can be a significant 
impact on the value of a 
plethora of device parameters 


such as depletions widths and/or 
FET threshold voltages if the 
assumption is violated. In the 
domain that th temperatur is 
too low the percentage 
ionization of dopant or dopants 











will be significantly less than 
100%. The value for the 
majority carrier concentration 


is depressed significantly below 
the design value. On the other 
hand if th temperatur is too 
high the thermal generation 
effect causes the majority 
carrier concentration to become 
excessively higher than the 
design value in what is called 
the intrinsic temperature 
regime. Th xhaustion regime 
lies between this two extremes, 
intrinsic and freeze-out. It is 
well known that for a multiple 











impurity dopant process this 
"simple" three-regime 
description can be inadequate. 
What is important to the 
designer is the plot of the 
majority carrier versus the 


temperature, or what is more 
commonly done, a plot of 
majority carrier concentration 
versus reciprocal of the 
temperature. 

To-date most "analytic" methods 
Lox determining dominant 
features in such plots make use 
of multiple ad-hoc arguments, 
which taken one at atime, 


applies in only two of the three 
regimes mentioned above. [5, 6, 
7] Numerical methods based on 
one-dimensional analysis are 
applicable and provide 
significant flexibility in terms 
of making predictions when there 
are multiple dopants and when 
taking into account 2°° order 
effects. [4] It can be argued 
that numerical methods are 
generally going to be applicable 
over a wider range of problems 











than what can be solved with 
exact analysis. Two and three 
dimensional numerical finite 
element methods provided in CAD 
based commercial packages have 
been widely mployed by the 
engineers in the semiconductor 
processing industry. The 1-D 
results both numerical and 
analytical are compared with 
predictions from such a package 
of tools from Silvaco Int. [1] 


See Figure (1). For details on 
numerical algorithm please see 
references [3] and [4]. 

2. Basic Analysis 
The physical basis [1] for 
analysis assumes that there is 


sufficient spatial uniformity in 
the doping profiles in order to 
apply the condition of local 
charge neutrality, 


P(T)= p,-n,- >\N, +>0N7 =0 (1) 








where T is the Kelvin 
temperature, p, and Ns are the 
equilibrium hole and electron 
carrier concentrations 
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‘ E ‘ = + 
respectively, while MN, and N, 
are the ionized acceptors and 
donors atom concentrations 


respectively. The summations run 





over the numbers of impurities. 
The approach presented is 
facilitated by defining the 
parameter Z where: 

Z = exp| (E,,-E,)/kT | (2) 


where k 
E, is th 
the intrinsic 


is Boltzman's constant, 
Fermi level and E,;; is 
Fermi level. In 
non-degenerate cases, i. wher 

the Boltzman approximation can 
be applied to simplify the exact 
expression for the carrier 
concentrations, the Z parameter 
is directly proportional to the 
equilibrium hole concentration. 
The constant of proportionality 
is th temperatur dependent 
intrinsic concentration. It 
turns out that Eq. (1) can be 
expressed as P(Z). Numerical 
methods based on successive 
substitution in solving for the 
zeros in expression (1) were 
found to be unreliable in 
producing a convergent solution. 
However the method proposed here 
is derived from the numerical 
scheme known as interval 
bisection [8] and it was found 
work well with variety of 
combinations of profiles 
other nonstandard conditions 









































Th ffective density of states, 
which is derived in most 
introductory semiconductor 
texts, is given by: 
2am kT)” 
Ng C= (3a) 
h 
« \3/2 
mM, 
N,=N.|) — (3b) 
mM, 
where m is the effective mass 
for th respectiv bands [5]. 
In order to facilitate casting 














Eq. (1) into a form dependent on 
th Z-parameter th well-known 
relations for intrinsic 
concentration and Fermi level 
are used. 
Be 
= 2kT 
n, = NN, e (4) 
EL.+k N 
E,, =——— + kT In| — (5) 
aN: 
Following the standard methods 


[5, 6, 7] of expressing carrier 





concentrations in terms of 
intrinsic parameters Eqs (4-5), 
it follows that: 
Nn; 
1,= > (6a) 
Z 
Pp, =n,L (6b) 
Making use of the Fermi 


probability distributions [5] to 
predict ionization levels of the 
donors and acceptors, 1.e., 


Nand N; : 
Eqs (6) 
dividing 
concentration 
condition on P 








and substitution of 
and Eq. (1) and then 

by the intrinsic 
[3] leads to a 
(Z, T): 














in which the over-bar notation 
indicates division by ni , and 
the temperature dependent 
constants, 

_ (Emu )/kT 
K, = 8 (8a) 

(E, —Ep;) (kT 

a ay Fi 
K, =e (8b) 
are needed to characterize 
partial ionization. Also M, and 
Mz, are the numbers of impurity 


components, 
Silicon) 


and g, 
and gg 


(4 typical for 
(2 typical for 
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Silicon) are the several 
occupation degeneracy, for 
acceptors and donors, 
respectively. The donor and 


acceptor concentrations with the 
over-bars are normalized by the 
intrinsic concentration. The 
energy levels for the impurities 
Bap Ba cate defined relativ to 
the valence band, i.e., 








E,=E,-E, (9) 


where E, is the standard cited 


value for donors measured with 
respect to the conduction band. 
This analysis also be the 
basis for an analytic 
description provided in the next 
section. 


will 





3. Analytic Solution, 
Single impurity 


Situation is described by the 
special case that M.,=1 and M,=0 








the summation subscript in Eq 
(7) can be dropped to lighten 
the notation. Straight forward 
algebra leads to a cubic 
equation in Z. 

Fete Gl ee a eG (10) 

K he K, 

The corresponding result for a 
single impurity N-type is 


exactly the same form, producing 
a cubic in Y=1/Z with revised 
coefficients obtained by letting 


the ‘a’ subscript (for acceptor) 
be replaced with the a6 led 
subscript (for donor). The 
solution can be defined in terms 
of coefficients for the 
reference cubic equation: 
3 2 = 
xX +ax° +a,x+a,=0 (1) 


To facilitate representation of 
the solution the following 
intermediate parameterization of 


the 
[9]. 


Q= (3a, -a?)/9 
R=(9a,a, — 27a, —2a;)/54 
D=Q°+R’ 


problem is commonly taken 


(12a) 
(12b) 


(12c) 


which is referred to 
“discriminant” for the 
cubic problem will dictate the 
type of solutions possible. As 
per the fundamental theorem of 
algebra there will be three 
roots. For D negative all roots 
are real while for D positive 
only one root is real while the 
other two are complex conjugate. 
If D is zero there will be 
repeated root. For the problem 
being evaluated it can be shown 
using symbolic mathematical 
methods that Q and D will always 
be negative. Furthermore it can 
then be shown that of the three 
roots only one is positive and 
therefore physically acceptable. 
That root will be predicted by 
the following recipe [9] 


where, D, 
as the 











L=*F=2 =0.00s(8/3)- 4 (13) 


where @=cos" (x/J-0") . Once 


Z is determined th 
carrier concentration, 
predicted from Eq (6b). 





majority 
Po, can be 





4. Silvaco Int. based 





Simulation 
Typically, semiconductor device 
problems of interest involve 
conditions such as external 
voltages and/or spatial non- 
uniformity in the semiconductor. 
Under these more general 
conditions, a solution can be 
obtained by solving Poisson’s 


equation. This is done with the 


Silvaco tool known as ATLAS [2]. 
It is possible to incorporate 
Fermi-Dirac statistics for 
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purposes of estimating 
incomplete ionization. 

Due to lack of space it is not 
possible to cover all the 
details in the algorithm 
involving Silvaco Tools. Here 
are some critical points. 
First, the Silvaco predictions 
are obtained from running the 
software over 2 dimensional grid 
of points covering a 2 
dimensional, uniform and 
equilibrium semiconductor. A 
center point is chosen in the 
grid and the majority carrier 


concentration is recorded. [This 
is done in an automated fashion 





for each temperature in a long 
list of temperatures covering 
the range of interest. The 
Silvaco based recordings 
described in the next section 


are obtained via this procedure. 
5. Examples 


Figure 2, which has on it 6 





plots of logl10 (Po) versus 
1000x(1/T) for various conditions 
(with and without band-gap 
narrowing (BGN) and from various 
sources (numerical, analytic, 
and Silvaco). The model for 
band-gap narrowing [6] used here 
TS 
al’ 
E,(1)=E£,, - (14) 


LB) 


The various curves are coded by 


number with the description in 
the legend. In Table L 
recommended values for Silicon 


including BGN parameters @ and 6 





are provided. Note, as 
expected, curves from the 
numerical and analytic methods 
are indistinguishable. 
Furthermore, as expected the 
impact of BGN is apparent in the 
curves at the higher 
temperatures. The difference 
between Silvaco predictions and 





the 1-dimensional modeling 













































































































































































efforts are slight but noted for 
furure anvest igat Pon. With and Without Bandgap Narrowing Plot 
Parameter Value 19.0000 
a 0.045 ev —1 Silvaco - w/o BGN 
Na 10° #/cm’ —2 Silvaco - w/ BGN 
Ego 1.llev 18.5000 |} ~~ 3 Analytic - w/o BGN 
(Silicon) —— 4 Analytic - w/ BGN 
ga 4 — 5 Numerical - w/o BGN 
BGN a [6] 0.00047 jsontd — 6 Numerical - w/ BGN 
BGN B [6] 636.0 2 
Table 1: Parameters used in 
Simulations $22000 
3&5 
6. Conclusions z ei 
4 17.0000 ¢ 
Both analytic one-dimensional Tr i 
models and numerical one- 
dimensional algorithms are 1&2 
essentially indistinguishable in 
terms of predictions. As ree 
expected, when BGN effects are 
included in the models the 
impact was observed at the eee 
higher temperatures. The 
difference between Silvaco 
predi ct i one bas ed two- ia Fe 2.500 4.500 6.500 8.500 10.500 
dimensional modeling and the ENE 
one-dimensional solutions were 
observed to be slight but worthy Figure 2: Composite Plot of 
of future investigation. Results 
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